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ABSTRACT 

Magnetohydrodynamic (MHD) turbulence is pervasive in astrophysical systems. Recent high- 
resolution numerical simulations suggest that the energy spectrum of strong incompressible MHD 

turbulence is E(k±) cx k^ 2 . So far, there has been no phenomenological theory that simultaneously 
explains this spectrum and satisfies the exact analytic relations for MHD turbulence due to Politano 
& Pouquet. Indeed, the Politano-Pouquet relations are often invoked to suggest that the spectrum of 
MHD turbulence instead has the Kolmogorov scaling —5/3. Using geometrical arguments and numer- 
ical tests, here we analyze this seeming contradiction and demonstrate that the —3/2 scaling and the 
Politano-Pouquet relations are reconciled by the phenomenon of scale-dependent dynamic alignment 
that was recently discovered in MHD turbulence. 

Subject headings: MHD — turbulence — solar wind — ISM: magnetic fields 



1. INTRODUCTION 

Magnetized plasma turbulence plays an essential 
role in many astrop hysical phenom e na, such as 
the solar wind (e.g. iGoldstein et al. I I1995D. inter- 
stellar scintillation (e.g. ILithwick fc Goldreichl 12001 . 
cosmic ray acceleration, propagati on and scattering 
in th e interste l lar m edium (e.g., iKulsrud fc Pearcd 
Il969t IWentzell I1974T) and thermal conduction in 
galaxy clusters (e.g.. iRechester fc Rosenbluth 119781 : 



IChandran fc Cowlevlll998HNaravan fc Medvedevll2001f ). 
The statistical properties of such turbulence can be in- 
ferred either indirectly from astronomical observations, 
such as scintillation of interstellar radio sources, or from 
in situ measurements, such as measurements of the mag- 
netic and velocity fields in the solar wind. Over a wide 
range of scales, turbulence in astrophysical plasmas can 
be modeled in the framework of incompressible magne- 
tohydrodynamics. When written in terms of the Elsasser 
variables the equations have the form 

(I +v A -V)z+(wV)z = -VP, (1) 

(i-v A .V)w+(z.V)w=-VP, (2) 

where the Elsasser variables are defined as z = v b 
and w = v + b, v is the fluctuating plasma velocity, b 
is the fluctuating magnetic field normalized by y/Aivpo, 
va = Bo/\/47rpo is the Alfven velocity corresponding 
to the uniform magnetic field Bo, the "total" pressure 
P = pj po + b 2 /2 includes the plasma pressure p and 
the magnetic pressure and po is the background plasma 
density that we assume to be constant. 

Current theoretical understanding of MHD turbu- 
lence largely relies on phen omenological m odels and 
numerical simulations (e.g., iBiskamd 120031 ). How- 
ever, there are certain exact results in the statisti- 
cal theory of turbulence that can be used to test 
these predictions. The results have attracted consider- 
able interest recently, due to a number of solar wind 



observation s where the exact relations were directly 
tested (e.g.JVasauez et al 1120071: iMacBride et 1171 120081; 
iMarino et al. 1 120081 : iPodesta et al. Il2009h . In what fol- 
lows we discuss the exact relations formulated for MHD 
turbulence, the so-called Politano-Pouquet relations, and 
we analyze to what extent they agree with the recent phe- 
nomenological and numerical predictions for the energy 
spectrum of strong anisotropic MHD turbulence. 

It is well known that the Kolmogorov theory for 
isotropic incompressible hydrodynamic turbulence yields 
the following exact relation for the third-order longitudi- 
nal structure fu nction of the velocity field in the inertial 
range (see, e.g.. lFrisch 111995t h 



(<5«i(r)) = --er. 



(3) 



Here 6vl(t) — [v(x + r) — v(x)] ■ r/r is the longitudinal 
component of the velocity difference between two points 
separated by the vector r and e is the rate of energy 
supply to the system at large scales. In a stationary 
state it coincides with the rate of energy cascade toward 
small dissipative scales and with the rate of energy dis- 
sipation. If one assumes that the fluctuations are not 
strong compared to the rms value of <Sv(r), and that 
Sv(r) ~ 5vl{t), one can dimensionally estimate from §3§ 
that (Sv 2 (r)) cx r 2 / 3 . The Fourier transform of the latter 
expression then leads to the Kolmogorov spectrum for 
the turbulent velocity field, E(k) cx fc~ 5 / 3 . 

Interestingly, analogous relation s hold for isotropic 
magnet ohydrodynamic turbulence. iPolitano fc PouauetJ 
(19983B) derived 



S. 



3L (r) ^ (Sz L (5w) 2 ) = --e« 



s. 



3L (r) = (6w L (5z) 2 ) = --e z r, 



(4) 
(5) 



where Szl and Swl are the longitudinal components of 
Sz and 6w, e w is the transfer rate of the w field and e z 



2 



is the transfer rate of the z field. If one now follows the 
analogy with the hydrodynamic case and assumes that all 
typical fluctuations scale in the same way (Szl oc Swl oc 
Sz oc Sw oc Sv oc 8b) one derives Sv r oc Sb r oc r 1 / 3 , which 
leads to the Kolmogorov scaling of the MHD turbulence 
spectrum. 

The results (|4I5|) can be extended to the case of 
anisotropic turbulence in the presence of a strong guid- 
ing field - a setting relevant for astrophysical tur- 
bulence where a large-scale field is always present, 
whet her due to external sources or large-scale eddies (see , 
e.g.. iMaron fc Goldreicbl 120011 : iMatthaeus et al. 1 119961: 
Milan o et al. II2001D . According to iPolitano fc Poucmetl 
(1998a|), the requirement of homogeneity allows one to 
derive the following differential relations in the inertial 
range of turbulence 



_d_ 

Q r l 



d_ 

dr 1 



{5z\5w) 2 ) -- 
{5w l {5z) 2 ) 



-Ae w , 



-4e z . 



(6) 
(7) 



Expressions (|4I5[) immediately follow if the correla- 
tion functions are assumed to be three-dimcnsionally 
isotropic. However, in the case with a strong guiding 
field the variations of the fluctuations in the field per- 
pendicular direction are much stronger than their field 
parallel variations. Hence the latter can be neglected in 
the inertial interval and the spatial derivatives in [7]) 
can be replaced by their field perpendicular parts 



d 



dr 



(5z\5w) 2 ) 



-4e w , 



— { S w *(Sz) 2 )=-4e* 



(8) 
(9) 



where r_i_ is a vector in the ficld-pcrpcndicular plane and 
the variations of the fields z and w are taken along r j_ , for 
example Sw = w(x + rj_) — w(x). Assuming statistical 
isotropy in the field-perpendicular plane one then derives, 
analogously to (|4I5[) 



(Sz L (Sw) 2 ) = 
(Swl(Sz) 2 ) 



-2e w r±, 
-2e z rj_. 



(10) 
(11) 



Here the longitudinal components of z and w are de- 
fined along the two-dimensional vector rj_, for example 
Szl = [z(x + rj_) — z(x)] • v±/r±. A more formal deriva- 
tion of expressions analo gous to (|10lll|) can be found in 
iPerez fc Boldvrevl (|2008l ). 

Arguments similar to those described for the hy- 
drodynamic case may lead one to conclude that the 
energy spectrum for MHD turbulence with a strong 

guide field is also E(k±) oc k~ 5 ^ 3 (see , e.g., the dis - 
cussion and references in iBiskampl 120031 : iVermal 120041) . 
Phenomenological arguments leading to such a spec- 
trum have attracted considerable attention ([Higdon 
fl98i IGoldreich fc Sridhari [19951 : IVermal l200i) . This 
however reveals a puzzling contradiction with recent 
high resolution numerical simulations of strongly mag- 
netized turbulence that instead suggest E(k±) oc 

k~ 3/2 (see[Maron fc Goldliicl[200ll: iM uller et al.l|200l; 
IMiiller fc Gr'apphl[200l iMason et al. I [2006L I2008T) . Al- 
though phenomenological argu ments explaining suc h a 
spectrum have been proposed (|Boldvrevll2005L 120061 ) so 



far it remained unclear whether they can be reconciled 
with the exact Politano-Pouquet relations. This appar- 
ent inconsistency motivated our interest in the problem. 

The goal of the present paper is to analyze whether 
the field-perpendicular energy spectrum E(k±) oc k ± 3 ^ 2 
is consistent with the exact Politano-Pouquet relations. 
We propose that the spectral exponent —3/2 does not, 
in fact, contradict the rigorous Politano-Pouquet result. 
Rather, their relation is manifested in the phenomenon 
of scale-dependent dynamic alignment that was recently 
discovered in driven MHD turbulence. 



DYNAMIC ALIGNMENT AND POLITANO-POUQUET 
RELATIONS 



Dynamic alignment is a known 
MHD turbulence (e. g.. iDobrowolnv et al.l 
2 iPoucmet et al. I 



Grappin et al. 



Politano et al. 



198 
I1989D . 



phenome n on of 
19801 : 



1986, 



1988; 



However, in previous studies 
it essentially meant that decaying MHD turbulence 
asymptotically reaches the so-called Alfvenic state 
where either v(x) = b(x) or v(x) = — b(x), depending 
on initial conditions. It has been realized recently 
that the effect is modified in randomly driven MHD 
turbulence, where the fluctuations S~v r and ±5h r tend to 
align their directions in such a way that the alignment 
angle becomes sca le-dependent (|Boldvrevl 120051 120061 : 
IMason et al~1l2006D . 

The essence of the phenomenon is that at each field- 
perpendicular scale r (~ l/fc±) in the inertial range, 
typical shear- Alfven velocity fluctuations (Sv r ) and mag- 
netic fluctuations (±5h r ) tend to align the directions of 
their polarizations in the field-perpendicular plane and 
the turbulent eddies become anisotropic in that plane. 
In such eddies the magnetic and velocity fluctuations 
change significantly in the direction almost perpendic- 
ular to the directions of the fluctuations themselves, <5v r 
and ±5h r , which reduces the strength of the nonlinear in- 
teraction in the MHD equations. A nu merical illustration 
of thi s phenomenon can be found in IPerez fc Boldvrevl 
(2009). The alignment and anisotropy are stronger for 
smaller scales, with the alignment angle decreasing with 
scale as 9 r oc r 1 / 4 . This leads to the velocity and mag- 
netic fluctuations 5v r ~ Sb r oc r 1 / 4 and hence the energy 
spectrum E(k±) oc fc ± 3 ^ 2 , discussed in the introduction. 
We will now show its effect on the Politano-Pouquet re- 
lations. 

There are two possibilities for the dynamic alignment: 
the velocity fluctuation <5v. r can be aligned either with 
Sh r (positive alignment) or with — Sh r (negative align- 
ment). This implies that the turbulent domain is frag- 
mented into regions of positive and negative alignment. 
If no overall alignment is present, the numbers of posi- 
tively and negatively aligned eddies are balanced on av- 
erage. In the case of nonzero overall alignment there is 
an imbalance of positive ly and negatively a l igned eddies. 

However, as found in IPerez fc Boldvrevl (|2009l ). there 
is no essential difference between overall balanced and 
imbalanced strong turbulence. Therein it is argued that 
strong MHD turbulence, whether overall balanced or not, 
has the characteristic property that at each scale it is 
locally imbalanced. Overall, it can be viewed as a super- 
position of positively and negatively aligned eddies. The 
scaling of the turbulent energy spectrum depends on the 
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Fig. 1. — Sketch of the velocity and the magnetic field fluctua- 
tions, <5v r and <5b r , aligned in the field-perpendicular plane within 
a small angle 9 r . 

way the alignment changes with scale, not on the amount 
of overall alignment. 

Let us now determine which of the two configurations 
(positive or negative alignment) provides the dominant 
contribution to the structure functions (J4j> and (J^. We 
note that, by definition, the amplitudes of <5v r and <5b r 
are of the order of their typical, rms values. In the case 
when S~v r is aligned with 6b r , we have Sw r > 5z r . There- 
fore, this configuration contributes more to structure 
function (j4]) than to structure function ((5|). Similarly, 
the configuration in which 5v r is aligned with —Sh ri and 
hence 5z r > Sw r , provides the dominant contribution to 
structure function (|5|). 

Without loss of generality we consider in detail only the 
structure function S™ L (r), defined in ([4]), and we concen- 
trate on the contribution provided by the configuration 
in which <5v r is aligned with Sh r inside a turbulent eddy. 
Figure [1] illustrates this case. The large-scale field Bo is 
in the z-direction and the vectors 8v r and Sh r are aligned 
within a small angle 9 r in the field-perpendicular plane, 
in the y-direction, say. Since the polarization of shear- 
Alfven waves are perpendicular to Bo and to their wave 
vector k, the wave vectors k are aligned in the x-direction 
here. Consequently, the variation of the fields is strongest 
in the x-direction and the dominant contribution to the 
structure function (j4|) comes from the situation in which 
the point-separation vector r lies in the x-direction. It 
follows from geometrical arguments that the longitudi- 
nal projection (i.e., x-componcnt) of Sz r , 5zl, is smaller 
than the typical value of Sw r by a factor of order 9 r , viz. 
Sw r ~ Sv r and Szl ~ 9 r 5v r . This introduces an extra 
factor of 8 r in Politano-Pouquet correlation function ((4]) , 
and one obtains 

(Sz L (Sw) 2 } ~ 9 r 6v 3 r . (12) 

As was de monstrated i n iBoldvrevi (|2005L [2006) and 
iMason et al. I (|2006ll2008t ). the scale-dependent dynamic 
alignment 9 r cx r 1 / 4 leads to the scaling of the fluctu- 
ating fields Sv r ~ Sb r oc r 1 / 4 , which explains the nu- 
merically observed field-perpendicular energy spectrum 

E(k±) oc k ± 3 ^ 2 . Quite remarkably, by substituting these 
scalings into expression (fT2"|) we also satisfy the scaling 
relation (UJ). Thus, the numerical findings are reconciled 
with the Politano-Pouquet relations if one invokes the 
phenomenon of scale-dependent dynamic alignment. 

In the next section we attempt to test relation (|12[) us- 
ing numerical simulations of strong incompressible MHD 
turbulence. In particular, we measure the following 
third-order structure functions 




<l6z L l (5w) z >/<l5 

<l6w L l (6z) ! >/<l<5zl 3 > 



Fig. 2. — The relative scaling of the structure functions S^ L (r) 
and S™(r). The ratio l|15| l is plotted versus scale r. The equivalent 
procedure with w and z interchanged yields a similar slope. The 
straight line has the slope 0.2. 

s 3 » = <|<H 3 )- (14) 

We use the absolute value of Szl in calculating (TT3]) to 
avoid cancellations and slow convergence caused by dif- 
ferent signs of 5zl . If our idea expressed by (TT2"]) is correct 

then the functions S™ L (r) and S 3 (r) should have essen- 
tially different scalings, and, as follows from the estimate 

Sv 3 , their ratio should cor- 



S 3 L ( r ) ~ and S 3 (r) 

respond to the alignment angle 9 r 
scale as 



and therefore should 



S 3 * L (r) = (\5z L \(5w)*), 



(13) 



^(r)/Sr(r)«r - 25 . (15) 

Note that the scaling of the alignment angle 9(r) cx r 25 
can be measured independentl y with the aid of second- 
order structure functions fe-g. . [Mason et al. Il2006l |2008; 
iBeresnvak fc Lazarian 1 12006 



3. NUMERICAL RESULTS 

We solve the incompressible MHD equations (JU [2|) us- 
ing standard pseudospectral methods. An external mag- 
netic field is applied in z-direction with strength Bo ~ 5 
measured in units of velocity. The periodic domain has 
a resolution of 512 3 mesh points and is elongated in the 
z-direction with aspect ratio 1:1 :Bq. An external force 
f(x, t), and small fluid viscosity v and resistivity rj are 
added to the equations. The external force is random 
and it drives the turbulence at large scales. The details 
of the numer i cal m ethod and set up can be found in 
IMason et al~l (|2008l ). 

The Reynolds number is defined as Re = U rma L/v, 
where L (~ 1) is the field-perpendicular box size, v is 
fluid viscosity and U rms (~ 1) is the rms value of velocity 
fluctuations. We restrict ourselves to the case in which 
the magnetic resistivity and fluid viscosity are the same, 
v = r], with Re w 2200. The system is evolved until a 
statistically steady state is reached, which is confirmed 
by observing the time evolution of the total energy of the 
fluctuations. The data set consists of 30 samples that 
cover approximately 6 large-scale eddy turnover times. 
To calculate the structure functions (Ti"3")) and (Ti"4"| we 
construct 8z(r) = z(x + r) — z(x) and <5w(r) = w(x + 
r) — w(x), where r is in a plane perpendicular to Bo- By 
definition, 8zl = 8z(r) ■ r/r and Swl — S~w(r) ■ r/r. The 
average is taken over different positions of the point x 
in that plane, over all such planes in the data cube, and 
then over all data cubes. 
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Fig. 3. — Alignment regions at different scales. Plotted is 
cos(<9) = (v ■ b)/|v||b| in a plane perpendicular to the guiding 
field Bo , where v and b are the velocity and magnetic field filtered 
at the scale kj. Left-right: kf = f,4, fO. 

The results we presen t here corresp o nd to one of five 
simulations described in lMason et al. I (|2008l ) (Case 2a). 
Each of those simulations differ by the large-scale driving 
mechanism, which is not expected to affect the turbulent 
dynamics in the inertial interval. The numerical calcu- 
lation of expression (fTS)) is shown in Figure [U For this 
case we find S™ L (r)/S™(r) cx r 2 . Repeating the calcu- 
lation for the other four cases yields noticeable scatter 
from case to case, with the slopes ranging from approx- 
imately 0.17 to 0.21. We discuss the results in the next 
section. 

4. DISCUSSION AND CONCLUSION 
We have proposed phenomenological arguments 

demonstrating that the energy spectrum E(k±) cx k^" 1 
for MHD turbulence is consistent with the exact analytic 
results that are known as the Politano-Pouquet relations. 
We argue that the two results are reconciled by the phe- 
nomenon of scale-dependent dynamic alignment. As a 
numerical illustration, we computed the structure func- 
tions (fj"3"l) an d (O- We found that their scalings are 
essentially different, which supports our arguments pre- 
sented in section [2l However, the agreement with the 
analytic prediction S™ L (r) / S™ (r) cx r 25 turns out to be 
only qualitative. Several reasons may explain the lack of 
good quantitative agreement. First, this may be due to 
extremely slow convergence of the statistics for the third- 
order structure functions in ([15]) . Indeed we noticed that 



in each run the convergence of the statistics for the mea- 
sured third-order structure functions is quite slow. The 
evolution of the slopes from snapshot to snapshot ap- 
parently has long-time variations. In spite of the large 
statistical ensemble accumulated in our runs, it is not 
enough for the precise measurement of the slope in (Tl~5|l . 
Second, the size of the Reynolds number in our simu- 
lations (which is limited by computational costs) may 
significantly impede computation of the slope. Third, 
there may be a systematic deviation of the slope mea- 
sured in (|15p due to intermittency corrections to the 
third-order structure function (j!4|L which are not cap- 
tured by our model. We plan to address these issues in 
future numerical work. 

Finally, in Figure [3] we illustrate the regions of align- 
ment at different scales. It appears that majority of the 
domain is covered with regions of high positive or neg- 
ative alignment, and that this structure is hierarchical 
- within a large-scale region of positive alignment, say, 
there exist smaller scale highly aligned and anti-aligned 
structures. This supports our assumption that even over- 
all balanced MHD turbulence is imbalanced locally, so 
that different spatial regions are responsible for differ- 
ent structure functions in the Politano-Pouquet relations. 
We note that the presence of correlated polarized regions 
in MHD turb ulence was also observed in some early sim- 
ulations (e.g. JMaron fc Goldreichll2001h . 
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